IHEP 2005-20 



February 7, 2008 
OEF 

Observation of the radiative kaon decay 

K~ — >• /i~7T°7Z/ 

O.G. Tchikilev, S.A. Akimenko, G.I. Britvich, K.V. Datsko, A.P. Filin, A.V. Inyakin, 
V.A. Khmelnikov, A.S. Konstantinov, I.Y. Korolkov, V.M. Leontiev, V.P. Novikov, V.F. Obraztsov, 
V.A. Polyakov, V.I. Romanovsky, V.I. Shelikhov, V.A. Uvarov, O.P. Yushchenko. 

Institute for High Energy Physics, Protvino, Russia 

V.N. Bolotov, S.V. Laptev, V.A. Duk, A.Yu. Polyarush. 

Institute for Nuclear Research Moscow, Russia 



Abstract 



Using data collected with the "ISTRA+" spectrometer during the 2001 run of the U-70 proton syn- 
chrotron in Protvino, we report the first observation of the radiative kaon decay K — > ^Tix°^v. We find 
BR(K M 3 7 , 5 <E; < 30 MeV) / BR(K m3 ) = 0.270 ± 0.029(stat) ± 0.026(syst)% and BR(K m37 , 30 <E; < 
60 MeV) /BR(K^ 3 ) = 0.0448 ± 0.0068(stat) ± 0.0099(syst)%. These ratios are consistent with the the- 
oretical predictions 0.21% and 0.047% respectively. The measured angular distribution asymmetry for 
the region 5 <E* < 30 MeV, A(cos6>* 7 ) = 0.093 ± 0.141, is two standard deviations away from the 
theoretical prediction of 0.354. The measured asymmetry in the T-odd variable £ = p 7 ■ (p^ x p^/m^- 
is -0.03 ±0.13. 

1 Introduction 

The study of the radiative kaon decays can give valuable information on the kaon struc- 
ture and allows for good test of theories describing hadron interactions and decays, like 
Chiral Perturbation Theory (ChPT). Until now the studies of the radiative K^-decays 
are restricted by the decay modes with electrons in the final state or by the studies of 
K L decays (HQ- Only one paper [3], dated by 1973, published the upper limit on the 
branching of K+ 37 -decay. 

The interest to the study of IQ3 7 decays is further enhanced by the theoretical pro- 
posals to search for effects of new physics using T-odd kinematical variable £ = 
p 1 ■ (pi x PTr)/Tn 3 K [01 13| . In the standard model the expected asymmetry for K~ 37 
decay 

a(c\ _ iv(e>o)-iv(£<o) 

^\<s) N(£>0)+N(£<0) 

is at the level 1.14 x 10 " 4 0], whereas in the extensions of the standard model it can 
achieve 2.6 x 1CT 4 Q. 

In this paper we present first observation of the radiative K~ 3 decay. The experi- 
mental setup and event selection are described in section 2. The results of the analysis 
are presented in section 3, where we, first, show the presence of the signal for the pho- 
ton energy in the kaon rest frame E* below 60 MeV, measure the branching ratio for the 
region 5 < E* < 30 MeV, measure the asymmetries in this region and finally measure 
the branching ratio for the region 30 < E* < 60 MeV. Our conclusions are given in the 
last section. 
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Figure 1: The elevation view of the ISTRA+ setup. Mi and M 2 are magnets, Q — Cerenkov counters, 
BPCj — beam proportional chambers, PQ — spectrometer proportional chambers, SPj — lead glass 
calorimeters, DQ — drift chambers, DTj — drift tubes, HC — a hadron calorimeter, Sj — trigger 
scintillation counters, MH — a scintillating hodoscope and MuH — a scintillating muon hodoscope. 



2 Experimental setup and event selection 

The experiment is performed at the IHEP 70 GeV proton synchrotron U-70. The IS- 
TRA+ spectrometer has been described in detail in recent papers on K e z [0 EJ, 
[[SI [9[ and 7T 7r°7r° decays [10]. Here we recall briefly the characteristics relevant to 
our analysis. The ISTRA+ setup is located in a negative unseparated secondary beam 
line 4 A of U-70. The beam momentum is ~ 25 GeV/c with Ap/p ~ 1.5%. The 
admixture of K~ in the beam is ~ 3%, the beam intensity is ~ 3 ■ 10 6 per 1.9 sec 
U-70 spill. A schematic view of the ISTRA+ setup is shown in Fig. 1. The beam par- 
ticles are deflected by the magnet Mi and are measured by four proportional chambers 
BPCi — BPC4 with 1 mm wire spacing, the kaon identification is done by three thresh- 
old Cerenkov counters Co — C2. The 9 meter long vacuum decay volume is surrounded 
by eight lead glass rings used to veto low energy photons. The 72-cell lead-glass 
calorimeter SP2 plays the same role. The decay products are deflected in the magnet 
M2 with 1 Tm field integral and are measured with 2 mm step proportional chambers 
PCi — PC3, with 1 cm cell drift chambers DCi — DC3 and, finally, with 2 cm diam- 
eter drift tubes DTi — DT4. The wide aperture threshold Cerenkov counters C3 ,Ca , 
filled with He, serve to trigger electrons and are not used in the present measurement. 
SPi is a 576-cell lead-glass calorimeter, followed by HC, a scintillator-iron sampling 
hadron calorimeter. MH is a 11x11 cell scintillating hodoscope, used to improve the 
time resolution of the tracking system, MuH is a 7x7 cell muon hodoscope. 

The trigger is provided by scintillation counters Si — S5, beam Cerenkov counters 
and by the analog sum of amplitudes from last dynodes of the SPi : T=Si • S2 • S3 • S4 • 

V V V 

Ci • C2 • C3 • S5 • S(SPi), here S4 is a scintillation counter with a hole to suppress the 
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beam halo, S 5 is a counter downstream of the setup at the beam focus, S(SPi) requires 
that the analog sum to be larger than the MIP signal. 

During the physics run in November-December 2001 350 million trigger events 
were collected with high beam intensity. This information is complemented by 150 M 
Monte Carlo (MC) events generated using Geant3 [fTT| for the dominant K decay 
modes, 100 M of them are the mixture of the dominant decay modes with the branch- 
ings exceeding 1 %, 30 M are the decays K~ — > fi~ir u (K^) and 20 M are the decays 
K — > 7r 7r°7r° (K^). Signal efficiency has been estimated using 5.7 M MC events of 
the radiative decay weighted with the matrix element, calculated in the leading 
approximation (up to terms of 0(p 4 )) of chiral perturbation theory [|4l [T2|| 

Some information on the data processing and reconstruction procedures is given in 
© IE1 HH El EL here we briefly mention the details relevant for present analysis. 

The muon identification (see fl8l|9]]) is based on the information from the SPi and the 
HC. The energy deposition in the SPi is required to be compatible with the MIP signal 
in order to suppress charged pions and electrons. The sum of the signals in the HC 
cells associated with charged track is required to be compatible with the MIP signal. 
The muon selection is further enhanced by the requirement that the ratio r% of the HC 
energy in last three layers to the total HC energy exceeds 5 %. The used cut values are 
the same as in 0. 

Events with one reconstructed charged track and three reconstructed showers in 
the calorimeter SPi are selected. We require the effective mass 01(77) to be within 
±40 MeV/c 2 from m^. In the following analysis the central ±20 MeV/c 2 band is used 
for signal search and the side bands 95 — 115 MeV/c 2 and 155 — 175 MeV/c 2 are used 
for background studies. We require also the reconstructed z-coordinate of the vertex to 
be below 1650 cm. 183672 events have been selected and written to miniDST's using 
the above cuts with relaxed cut on to be above 1 %. 

3 Evidence for signal and measurements of the branching ratios 

A set of cuts is developed to suppress backgrounds to the K M 3 7 decay and/or to do data 
cleaning: 

0) We select events with good charged track having two reconstucted (x — z and 
y — z) projections and the number of hits in the MH below 5. We require also that 
the missing mass squared to the (^ 7r°7) system abs(m 2 (/i~7r°7)) < 0.05 (GeV/c 2 ) 2 . 
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50804 events have survived these cuts. 

1) Events with the reconstructed vertex inside the interval 400 < z < 1600 cm are 
selected. 

2) The measured missing energy E mis = Eb eam — — E n o — E 1 is required to be 
above zero. 

3) We require the effective mass M(77) to be within ±20 MeV/c 2 from nv. 

4) We require also that the missing mass squared to the 7r 7r° system is below 
0.025 (GeV/c 2 ) 2 (m 2 (7r-7r°) < 0.025). 

5) The events with missing momentum pointing to the SPi working aperture are 
selected in order to suppress possible 7r _ 7r°7 background ( 6 < r < 60 cm, here r is 
the distance between the impact point of the missing momentum and the SPi center in 
the SPi transverse plane). 

6) We require the photon energy E* in the kaon rest frame to be below 60 MeV. 
The remaining decays are suppressed by requirements: 

7) cos(0) > —0.96 , where 6 is the angle between 7r and tt° in the kaon rest frame; 

8) (p < 3.0, where (p is the angle between 7r and tt° in the laboratory frame in the 
plane perpendicular to the beam momentum. 

9) We require also the absence of the signal above the threshold in the calorimeter 
SP 2 . 

We look for a signal in the distributions over the effective mass M(/i~7r°7z/), where 
v four-momentum is calculated using the measured missing momentum and assuming 
m v = 0, and in the distributions of the missing mass squared to the (/i~7r°7) -system, 
m 2 ([i~7r 0/ -f). Effective mass spectra for cut levels 1, 4, 6 and 9 are shown in Fig. 2. 
These spectra show the evidence for peak at m K after the cut on the photon energy in 
the kaon rest frame. 

3.1 The region below 30 MeV. 

We have found that the signal is clearly seen for E* < 30 MeV and the background in 
this region is dominated by decays (with an accidental extra photon) and K^-3 de- 
cays. The main MC sample of 100 M events (with the natural mixture of the dominant 
decay modes) has been found to be insufficient for estimates of the background shapes, 
therefore specialized MC samples of 20 M and 30 M K^ 3 events have been used. 
The background has been divided into three contributions: 
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Figure 2: Effective mass M(/x 7r°7z/) for the cut levels 1,4,6 and 9 respectively. 
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1) Non-7r° contribution has been estimated using tails of the M(77) distribution for 
real data, see Fig. 3 . 

2) contribution has been approximated by the form given by specialized MC 
sample, its normalization has been fixed using the observed signal in the m 2 (7r 7r°) 
distribution for selected events. 

3) K^3 contribution has been approximated by the form from specialized MC sam- 
ple, its normalization has been kept free. 

The shapes for all three background contributions have been found using the his- 
togram smoothing by the HQUAD routine from the HBOOK package IfT^ll . The signal 
has been parametrized by the sum of two Gaussians with widths and relative fractions 
fixed at the values given by the signal MC sample. 

Results of the fits are illustrated respectively in Fig. 4 and Fig. 5 for the distributions 
over M(httjv) and m 2 (/i7T7). First parameter here (and in the following) is the number 
of observed events, second parameter is the position of the peak , and three last param- 
eters are the respective normalization factors of the K^, and non-7r° contributions. 

The number of observed events is equal to 383.7 ± 40.9 in Fig. 4 and to 412.9 ± 
36.2 in Fig. 5. The difference between these two values (29.2) is our estimate of the 
systematics caused by the imprecise knowledge of the backgrounds. The results of fits 
with the polynomial parametrization of the background lie also within this uncertainty. 

The K^3 decay has been used for the normalization. The number of K M 3 events 
in the region 400 < z < 1600 cm, corrected for the geometrical acceptance and the 
measurement efficiency, has been found using the cuts described in [9]. It is equal to 
N(K^) = 5536000. Independent normalization of using decays gives the branch- 
ing ratio lower by 6.2%. The origin of this difference is explained mainly by the trigger 
bias. This difference has been taken into account in our final estimates of the systematic 
uncertainties. 

The signal efficiency has been found from the signal MC weighted with the matrix 
element, calculated within 0(p A ) ChPT approximation. It is equal to 2.6 %. The signal 
efficiency has been calculated using the cut E* > 5 MeV, the cut value is our detection 
threshold explained by the beam momentum and the threshold in the energy of the SPi 
showers equal to « 0.5 GeV. 

The measured branching ratio is equal to BR = (8.82 ± 0.94(stat) ± 0.86(syst)) x 
10 5 . This should be compared with theoretical prediction of 6.86x10 5 . The ra- 
tio R=BR(K M 3 7 )/BR(K / , 3 ) is equal to R = (2.70 ± 0.29(stat) ± 0.26(syst)) x 10~ 3 . 
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Figure 3: Effective mass M(77) for events with E* < 30 MeV at the cut level 9. Solid line shows 
parametrization by the sum of the Gaussian and polynomial background. 
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Figure 4: Effective mass M(/i - 7r°7z/) for events with 5 <E* < 30 MeV. Cross-hatched area shows the 
summary background, left(red) bump shows K„-3 contribution, right(blue) bump shows contribution. 
First parameter here and in the following figures is the number of signal events, second parameter is the 
position of the peak, third parameter is the normalization factor for contribution, fourth parameter is 
the normalization factor for K M 3 contribution and fifth parameter is the normalization factor for non-7r° 
contribution. 
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Figure 5: Missing mass squared m 2 (/i _ 7r°7) for events with 5 <E* < 30 MeV. Cross-hatched area 
shows the summary background, right(red) bump shows K^ 3 contribution, left(blue) bump shows 
contribution. 
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This should be compared with theoretical prediction of 2.1 x 10 3 . In the transforma- 
tions from the ratio R to the branching ratio we use the branching ratio BR^K^) = 
3.27% El. 



3.2 Asymmetries for the region 5 <E* < 30 MeV. 

For this region we have measured the asymmetry of photon emission towards muon 
direction in the kaon rest frame A(cos #* 7 ) 

AfWjfl* ^ - iV(co S ^ 7 >0)^iV(cos^ 7 <0) 
/^UJb u^) AT(cos0* 7 >O)+iV(cos0* 7 <O) 

and the asymmetry in £, A(£). The effective mass spectra are shown in Fig. 6 for 
positive and negative cos (9* 7 separately and in Fig. 7 for positive and negative £. 

The asymmetry in cos #* 7 is equal to 0.09 ± 0.14, this value is below the theoretical 
expectation, equal to 0.35, by 2 standard deviations. 

The asymmetry in £ is equal to —0.03 ± 0.13. Of course our statistics is insufficient 
to test the theoretical predictions 015]]. 



3.3 The region 30 <E* < 60 MeV. 

For this region we see strong background and residual K n2 background. These 
backgrounds have been suppressed by additional cut 0.1< p*(tt~) < 0.185 MeV/c. 
The parametrization is illustrated in Fig. 8 for the effective mass M(/z _ 7r°7^) spectrum. 
The signal efficiency was found to be 6.30 %. The number of observed events is equal 
to 152.7 ± 23.06(stat) ± 32.2(syst). The systematics in the number of observed events 
has been calculated in the same way as in the Section 3.1. The branching ratio BR = 
(1.46 ± 0.22(stat) ± 0.32(syst)) x 10 5 , is compatible with the theoretical expectation 
1.53 x 10~ 5 . The ratio R = (4.48 ± 0.68(stat) ± 0.99(syst)) x 10^ 4 , is compatible with 
the expectation 4.67 x 10 4 . 
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Figure 6: Effective mass M(/i 7T°7z/) spectra for: a) cos( 
background notations are the same as in Fig. 4. 
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Figure 7: Effective mass M( / u 7r°7z/) spectra for: a) £ > and b) £ < 0. The background notations are 
the same as in Fig4. 
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Figure 8: Effective mass M( / u 7r°7z/) spectrum for 30 <E* < 60 MeV. Cross-hatched area shows the 
background. Solid areas, left(red) and rigth(blue), show and K /t3 contributions respectively. 
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4 Conclusions 



Our conclusions are the following. 

• First observation of the radiative kaon decay K~ 37 is presented. 

• The measured ratio R = BR(K Al37 ) /BRCK^) for the region 5 <E* < 30 MeV is 
equal to 0.270±0.029(stat) ± 0.026(syst)%. This is consistent with theoretical predic- 
tion equal to 0.21 %. 

• The measured ratio R for the region 30 <E* < 60 MeV is equal to (4.48 ± 
0.68(stat) ± 0.99(syst)) x 10 4 , this value is compatible with theoretical prediction 
equal to 4.67 x 10 4 . 

• The measured asymmetry in the T-odd variable £ for the region 5 <E* < 30 MeV 
is equal to -0.03 ±0.13. 

• The measured asymmetry in the cos 6** 7 is equal to 0.093 ± 0.141, this value is two 
standard deviations away from the theoretical prediction equal to 0.354. 

The work is supported in part by the RFBR grant N03-02-16330(IHEP group) and 
RFBR grant N03-0216135(INR group). We are indebted to V. Braguta for giving us 
routine for 0(j? 4 ) ChPT matrix element calculations. 
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